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Abstract. We develop the gradient expansion formalism for shift-symmetric Galileon-type 
actions. We focus on backgrounds that undergo Galileon inflation, work in the synchronous 
gauge, and obtain a general solution up to second order without imposing extra conditions at 
first order. The solution simplifies during the late stages of Galileon inflation. We also define 
a curvature perturbation conserved up to first order. 

Keywords: Modified Gravity, Cosmological Perturbation Theory, Inflation, Non-Gaussianity. 



Contents 

1 Introduction 1 

2 Inflation with a shift-symmetric Galileon 3 

3 3 + 1 decomposition 3 

4 Gradient expansion: order analysis 5 

5 General solution 7 

5.1 The O(e ) order 7 

5.2 The 0(e) order 7 

5.3 The C(e 2 ) order 9 

5.4 Summary 11 

6 Late time of Galileon inflation 12 

7 Conclusion 14 

A Dependence of equations of motion in general covariant scalar-tensor the- 
ory 15 

B Gradient expansion of the equations of motion in synchronous gauge 15 



1 Introduction 

Inflation is a powerful paradigm that addresses various fine-tuning problems in the early 
Universe and accounts for the nearly scale invariant primordial perturbations that are needed 
for structure formation. These primordial perturbations leave an imprint on the cosmic 
microwave background (CMB). Single field slow-roll inflation models (with a canonical kinetic 
term) generically predict a Gaussian spectrum of primordial perturbations [1]. However, 
despite the success of the inflationary paradigm, its theoretical underpinning is still a matter 
of debate. Hence, it is not clear why one should remain within the framework of slow- 
roll inflation or even single-field inflation. Substantial non-Gaussianity can be generated in 
inflation models with multiple scalar fields or with non-canonical kinetic terms. Furthermore, 
if the slow-roll condition is temporarily violated, large non-Gaussianity can be generated even 
in a single field model [2]. On the observational side, results from WMAP are consistent with 
a Gaussian spectrum of primordial perturbations [3] and the recently released results from 
PLANCK are already leading to tighter constraints of non-Gaussianity [4]. With the prospect 
of probing the inflation scenario much deeper, non-Gaussianity in different inflation models 
has been extensively investigated and classified during the last decade (for a review, see [5, 6]). 

To tackle non-Gaussianity from inflation models, traditional linear perturbation theory 
is inadequate. A natural approach is to go beyond the linear order and work with second 
order cosmological perturbation theory [1, 7-12]. While this approach usually applies to 
primordial perturbations up to the horizon exit, an alternative approach naturally tackles 
the superhorizon perturbations - gradient expansion [13-25]. In gradient expansion, physical 
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quantities are expanded in terms of their inverse wavelengths, as compared to a pivotal 
length scale (e ~ L p / L p h ys ) 1 so every spatial derivative adds one perturbative order, di ~ 
e, hence th e name. This is different from usual cosmological perturbation theory where 
the expansion is in terms of perturbative field amplitudes. In the context of cosmology, 
particularly in the inflationary epoch when physical modes are stretched well outside the 
quasi-constant Hubble horizon, the Hubble length can be naturally chosen as the pivotal 
length scale. Therefore, this approach can be used to evaluate and evolve non-Gaussiantites 
at super horizon scales, complementary to usual non-linear perturbation theory. The leading 
order gradient expansion is often called the separate universe approach [26] or 5N formalism 
[27, 28], which is sufficient for many purposes. However, the next-to-leading order gradient 
corrections can be as important, for example, in some multi-field models or when the slow-roll 
condition is violated [25, 29]. A beyond-5iV formalism scheme has recently been proposed 
[25]. 

The Galileon is a scalar field with a galilean(-like) symmetry (f> — >■ <j) + b^x^ + c {b^,c 
being constant) around flat space, and was originally introduced as an effective, infrared 
gravitational modification which can lead to self-accellerating solutions [30]. The Galileon 
Lagrangian contains higher order derivatives but nevertheless leads to second order equations 
of motion, thus avoiding Ostrogranski ghosts. 

Coupling the Galileon covariantly to gravity and insisting on the requirement that the 
scalar and the metric satisfy second order equations forces one has to abandon galilean sym- 
metry [31]. Ordinary shift symmetry <fi — > <p + c can be retained [32] or abandoned as well. 
In the latter case, covariant actions for generalized Galileons coupled to gravity have been 
constructed in Refs. [33, 34]. The 4-dimensional version of the action given in Ref. [34] has 
been shown [35] to be equivalent to the most general action for a scalar coupled to gravity 
that leads to second order equations of motion, given by Horndeski in the 1970s [36]. The 
Galileon model has also been generalized in various other directions (see e.g. [37-43] and 
references therein) . 

The self-accelerating solutions of Galileon models have been the basis of inflationary 
scenarios [32, 35, 44-50]. There are some known novel features in Galileon inflation (or 
G-inflation [32]): The null energy condition can be drastically violated without developing 
instabilities [44]; A large tensor-to-scalar power spectrum ratio is allowed [32]; There are new 
shapes of the three-point function and potentially large four-point function [45]. 

In this paper, we develop the super horizon gradient expansion formalism for Galileon 
inflation, up to second order in gradient expansion. We focus the subclass of actions for which 
the scalar enjoys shift symmetry, as they are closer to the original idea of the Galileon and sig- 
nificantly simpler. Additionally, our goal is to explore the phenomenology associated with the 
non-linear derivative interactions of the scalar. Abandoning shift symmetry leads, amongst 
other terms, to allowing a potential for the scalar, which can lead to similar phenomenology 
and, therefore, obscures the role of the Galileon-type terms. We derive the general solution 
for an Friedman-Lemaitre-Robertson-Walker (FLRW) background, identify the degrees of 
freedom in the model and define a curvature perturbation conserved up to 0(e). We also 
consider how our results simplify in the limit where the background becomes de Sitter space- 
time. Given that the latter is an attractor of G alileon inflation, this approximation provides 
a good description at least for the later stages of inflationary expansion. 

During the preparation of this manuscript ref. [24] appeared, which also develops the 
superhorizon gradient expansion formalism for G-inflation without assuming shift symmetry. 
However, there are major differences: On the technical side, we work in the synchronous 
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gauge, while ref. [24] prefers the uniform expansion gauge; On the more substantial side in 
ref. [24] it is assumed that dthij(t,x) (see eq. 4.3) is 0(e 2 ), while we do not impose such 
condition. In this respect our results are more general. 

The paper is organized as follows. In Section 2, we specify the model to investigate 
and write down the equations of motion and in Section 3 we perform 3 + 1 decomposition of 
the equations of motion. In Section 4, we establish the gradient expansion orders of relevant 
quantities. The equations of motion are solved up to 0(e 2 ) in Section 5, and the general 
solution is summarized in Section 5.4. In Section 6, we simplify the general solution in the 
de Sitter limit. Note that in Section 5, Section 6 and Appendix B, we mostly suppress the 
background quantities' order indication to simplify the equations. Section 7 contains our 
conclusions. 



2 Inflation with a shift-symmetric Galileon 

Following the arguments in the Introduction, we consider the action 



/ 



M 2 \ 
d 4 x^j \-f-R + K(X)-G(X)D<f> , (2.1) 



where M p i is the Planck mass, R is the Ricci scalar and K(X) and G(X) are general functions 
of X = —d fJi (pd fl (j)/2. Note that for G = we recover k-inflation [51], and the model with 
K = X, G = aX (a =const.) is the cubic covariant Galileon [31], which in turn is linked 
to the DGP model [52]. The equations of motion for the metric, to which we will refer as 
Einstein equations, are given by 

M 2 pl G, u = T^ (2.2) 

with 

1% = (K x - GxU^d^d^ - 2d [lx Gd v) (t> + g^(K + d a G8^), (2.3) 

where a subscript X denotes partial differentiation with respect to X. Note that the energy 
momentum tensor takes the form of an imperfect fluid, thus this model does not fall under 
the existing formalism for a perfect fluid [20]. Thanks to the shift symmetry, the equation of 
motion for <j) can be given in terms of the current 

= (K x - G x a^(t> - G x d^X, (2.4) 



as 

V^J" = 0. (2.5) 

It is worth mentioning that the scalar equation of motion is implied by the Einstein equations, 
i.e., once the Einstein equations are satisfied, the scalar equation of motion is automatically 
satisfied. In fact, this applies to any covariant scalar-tensor system, as we show in Appendix A. 



3 3 + 1 decomposition 

Now, we perform 3+1 decomposition of the equations of motion. First, we decompose the 
metric according to the Arnowitt-Deser-Misner (ADM) prescription 

ds 2 = -N 2 dt 2 + (dx* + N*dt) (dx j + N j dt) . (3.1) 
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To reduce redundant gauge degrees of freedom and simplify equations, we make use of a gauge 
condition: 

N = 1, N* = 0, (3.2) 

which implies 

9tt = -l, gti = 0, = 7ij, (3.3) 
5 « = -l, / = 0, ^'=7°'- (3.4) 

Here latin indices (except for t) run from 1 to 3 and they are raised and lowered with ^ and 
7ij respectively. This is called synchronous gauge, where the proper time distance between 
two neighboring hypersurfaces along the normal vector coincides with the coordinate time 
distance (JV = 1, proper time slicing) and the spatial coordinates are such that clocks are 
synchronized between different hypersurfaces (N l = 0). In synchronous gauge, equations can 
be very much simplified. Note, however, that there is residual gauge freedom, which will be 
discussed in section 5.4. 

Now, the spatial 7^ can be considered as a fundamental dynamical variable. Another 
fundamental variable after the 3 + 1 decomposition is the extrinsic curvature, which in syn- 
chronous gauge is simply 

and its trace is defined as K, = j lJ fCij. It is useful to decompose the spatial metric and the 
extrinsic curvature as 

lij = a 2 (t)e 2 ^h ij (t,x), (3.6) 
q = iAC(t,x)dj. + 4(t,x), (3.7) 

where a(t) is the scale factor of the fiducial FLRW background, C(^ x ) is related to the 
curvature perturbation, hij(t,x) is defined to have a unit determinant det[/ijj] = 1, and A 1 - is 
the traceless part of /Q. These definitions lead to the following relations 



3 

' a 



K. 



+ c 



(3.8) 



a 

hij = ~2h ik Al (3.9) 

To decompose the equations of motion, we first make use of some well-known results 
which do not make reference to any specific gauge. Using the ADM variables, the unit 
normal 1-form and vector can be written respectively as = (-N, 0, 0, 0) and = 
(l/N, -N l /N, -N 2 /N, -N 3 /N). Making use of the Gauss-Codazzi relations (see e.g. [53]), 
we can write the Ricci tensor and Ricci scalar respectively as 

iV = n„n v (JjC m lC + ^T> X D X N - - 2n (fl D v) K + 2n (M |D (T /C£ ) 

- " -^ D ^ + [3li V + K*W " 2/C£/C„ ff) (3.10) 

R = Wr + K 2 + - ^£ m /C - ~D°D a N, (3.11) 



-A- 



where = Nn^, C m is the Lie derivative along m^, and is the covariant derivative, PI J? 
the Ricci scalar, and ^R^ v the Ricci tensor of the spacelike hypersurfaces. The Laplacian is 
decomposed as 



+ Krfdvcj) + D CT In Nd a <p + D a D CT 0. 


(3.12) 


N l = 0), the Einstein equations are . 


greatly simplified: 


MpiG tt = T tt , 


(3.13) 


MpiGti = Tf i: 


(3.14) 


M^G) = T*\), 


(3.15) 



with 



G tt = l -^R + K 2 -V ] lCl), (3.16) 

G ti = -D fc /C, fc + Di/C, (3.17) 

G) = ^G) - K) + KK) - \b) (-21C + K 2 + /CfJCjQ , (3.18) 

2$ = K x <j> 2 -K- GxOcjxj) 2 - G X 4>X - G x d k Xd k 4) : (3.19) 

t£ = Kxidit - G x n^di4> - G x xd i( f> - G x diXj>, (3.20) 

T*|j = Kxd^djcj) - G x D0</>dj<f> - Gx&Xdrf - Gx&WjX 

+ (G x d k <\>d k X + K- G X X$) 5), (3.21) 



where □</> = —0 + /C0 + D CT D CT 0. The scalar equation of motion is given by 

<V M + ify ln 7 J^ = J* + (3# + 3C) J* + 3ftC^ + ft J i = 0, (3.22) 
where 7 = det[ 7 y] and H = a/a is the usual Hubble parameter. 

4 Gradient expansion: order analysis 

In standard cosmological perturbation theory one expands perturbatively in the field ampli- 
tudes. To tackle non-Gaussianities in inflation models, second order perturbation theory is 
often used within the Hubble horizon. However, for physics processes at superhorizon scales 
one usually resorts to the gradient expansion technique. Note that the separate universe 
approach or the 5N formalism is simply the leading order gradient expansion [54]. Assuming 
the characteristic spatial length is L p h ys , the dimensionless perturbative expansion parameter 
is 

TT-l 

e~- <1. (4.1) 

l-"phys 

This means in particular that every spatial partial derivative carries an order of e 

ft ~ 0(e), (4.2) 

while the time derivative is considered 0(e ). The superhorizon gradient expansion is com- 
plimentary to the usual non-linear cosmological perturbation analysis and may capture fully 
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nonlinear (in terms of the field amplitudes) physics at superhorizon scales, while the equations 
are still tractable due to the perturbative approach. 

To perform the superhorizon perturbation analysis, we first need to deduce the starting 
orders for various quantities of interest. First, note that the equations of motion at 0(e°) 
should simply determine the evolution of the scale factor a(t) and the scalar, as the spacetime 
is supposed to be described by an FLRW line element. Given the definition (3.6), one can 
infer that hij should start with hfj ' (x); otherwise the 0(e°) equation would pick up terms 

involving dth^ (i,x), which is non-FLRW. From the scalar's equation of motion eq. (3.22), 
we can infer that at 0(e°) the scalar field should be spatially homogeneous, meaning that <fi 
starts with <fi(°'(t). Unlike previous work on this subject (see e.g. [24]), we do not impose any 
conditions on the higher orders of these quantities. Therefore, we have 

starting order of hij = 0(e), (4.3) 
starting order of 3$ = 0(e 2 ). (4.4) 

Expanding eq. (3.9) perturbatively (for n > 1) 

1 u 

and making use of eq. (4.3), we can infer that 

starting order of Aj = 0(e). (4.6) 
Expanding eq. (3.14), we infer that 

diK^ = 0. (4.7) 

Therefore, is a function of t. From the definition (3.8), and given that one can always 
redefine the scalar factor a(t) to absorb C it follows that 

starting order of ( = 0(e), (4.8) 

/C<°) = -3- = -3H(t), (4.9) 
a 

K.W = -3C (n) , n>l, (4.10) 



where H(t) is the usual Hubble parameter. 
Using eq. (4.4), we may expand X as 

X = X (0) (t,x) +X«(i,x)e + X( 2 )(i,x)e 2 + 0(e 3 ), (4.11) 

where 

1 n 

x (n) = 0(p)0(«-f) + Q{e A ). (4.12) 

p=0 

We also need to perturbatively expand functions of X, such as K(X). To this end, we should 
consider X = X(e) according to eq. (4.11) and Taylor-expand, for example, K(X(e)) around 
e = as 



K(X) = K(X<®) + K x (X^)X^e + \ [k xx {X^)X^ 2 + 2K X (X^)X^ 



e 2 + 0(e 3 ). 
(4.13) 
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In summary, the various quantities of interest, to be determined in the next section, are 
expanded as follows: 



c 


= C (0) (x) + C (1) (t, x)e + C (2) (t, x)e 2 + 0(e 3 ) , 


(4.14) 




= <A (0) (t) + (1) (t,x)e + </> (2) (t,x)e 2 + 0(e 3 ) , 


(4.15) 


4' 


= A^)(t, x)e + A< 2 >}(i, x)e 2 + 0(e 3 ), 


(4.16) 


h- 


= fc{°>(x) + h${t,x)e + hf(t,^)e 2 + 0(e 3 ) , 


(4.17) 


/q 


= -H(t)8) + KP>% x)e + x)e 2 + 0(e 3 ), 


(4.18) 


K 


= -3#(i) + C (1) (t, x)e + C (2) (t, x)e 2 + 0(e 3 ) . 


(4.19) 



5 General solution 

Now, we solve the equations of motion perturbatively, up to the 0(e 2 ) order, to obtain the 
general solutions. These solutions will be parametrized by a few unspecified spatial functions, 
which describe the physical degrees of freedom (modulo residual gauge freedom) that may 
evolve as the Universe expands. The gradient expansions of Einstein's equations up to order 
0(e 2 ) are listed in Appendix B. 

For the rest of this paper, to simplify the equations, we will mostly suppress the back- 
ground quantities' order indication For example, <^°) is written as (j) if there is no confu- 
sion. 

5.1 The 0(e°) order 

For the O(e ) order, all spatial derivatives are absent. As desired, the equations of motion 
reduced to the conventional background FLRW case: 

3M 2 t H 2 = -K + 2K X X + 6H<j)G x X, (5.1) 

-M 2 pl (2H + 3tf 2 ) = K — 2G X X4>, (5.2) 

+ 3H J*(°) = 0, (5.3) 

where = K x <\> + 6HG X X. Note that only two of the three equations are independent. 

5.2 The 0(e) order 

The tt component of Einstein's equation is given by 

^K x + <PK XX X + 9HG X X + 6HG xx X 2 ^j = (3M*H - 30G x a) C (1) , (5.4) 
which can be re-written as 

^=A°C {1 \ (5.5) 

where 

6M 2 H - 6<j>G x X 

A°(t) = -. = p - — . (5.6) 

<j)K x + 2<j)K xx X + 18HG X X + 12HG XX X 2 

The ij component of Einstein's equation naturally splits into a trace part and a traceless 
part. The trace part gives rise to another relation between </>W and £W: 

- 2G X X^ + Uk x - 2X(G X + G XX X)) jF* + 2M 2 (c (1) + 3#C (1) ) = 0. (5.7) 
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Combining eq. (5.5) and eq. (5.7), we get, after integration, 



df 



(5.8) 



where 



a(t) = exp( /" dt'H°{t')), (5.9) 

*"<«> - " m^tbuv -■ (5 - 10) 

and Cj 1} (x) is an unspecified spatial function from the first integration. There would be 
another unspecified spatial function from the second integration (C^ (x)), which has been 

absorbed into £(°)(x). We will see in the next section that Ti° approaches the Hubble constant 
H for near de Sitter expansion, in which case Q 1 ^ scales as l/a 3 (i). Now, is given by 



(1) (^) = C«(x)/ a ^ + C«(x), (5.11) 



a(t> 

where C, (x) is an integration spatial function. The traceless part of Einstein's equation's 
ij component is simply 

+ 3HA^i = 0, (5.12) 

whose solution is 

#H( t ,x) = ^, (5.13) 

where the unspecified spatial function C^* (x) is symmetric and traceless. From eq. (3.9), 
we have 

ft df 



Mj ) (^) = -2^)(x)C^(x)/-| : 



3 ' 



(5.14) 



where the would-be integration spatial function C^'^fx) has been absorbed into /i^(x). As 
expected, the scalar equation of motion is identically the solution obtained above. 

Before moving on to solve higher order equations, we note that defining a curvature per- 
turbation that is conserved in time is trivial in our formalism. By virtue of the tt component 
of Einstein's equation (5.5), one can define a conserved curvature perturbation at 0(e) 

* W = < W " W)- (5-15) 

As we will see in Section 6, because of the shift symmetry, de Sitter expansion is an attractor 
of the system. For quasi-de Sitter expansion, i.e., for the late time of Galileon inflation, 
A — constant and we can write 7iP~' as 

ftM~C (1) -^ (1) - (5- 16 ) 
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^(2) _ ^2) = C ( ^ _ ^ m R (2) + ^ A {l)k A (l)^ (5 17) 



5.3 The 0(e 2 ) order 

The tt component of the Einstein equation gives 

°3[3]r(2) + ^ 

2 2 

where I 3 li?( 2 ) is the 3D Ricci scalar for the O(e ) order metric 7^ = a(i) 2 e 2 ^ ( °^ x - ) /i^(x) and 
C°, C®, C\ and Cg are again background quantities, defined respectively as 

C°(t) = C?^ ) 2 + dU° - 3C§, (5.18) 

C°m = ^ Kx + 4J ^ xX + 2K ^xX 2 + 9ff</>G x + 21ff0G xx X + 6H<pG xxx X 2 
1 c^Kx + 2<j)K xx X + 18HG X X + 12HG XX X 2 

C 2 °(t) = , 18GxX + 12Gxx* 2 

+ 2(f)K xx X + 18HG X X + 12i/G X x^ 2 
M 2 

CS(t) = ; v - . (5.21) 

(j)K x + 24>K XX X + 18HG X X + 12tf G X x^ 2 

Integrating this equation, we get the solution of (jffl in terms of ^ 2 ^: 

^ = / W)C^,x) - (<f w)7 ^ + ^ / ^ 

2 J (5 - 22) 

where an integration spatial function has been absorbed into Cx (x), and [ 3 ]i?( 2 )(x) (the 
Ricci scalar of the metric e 2 ^°^ x )/iy (x)) is related to ^he Rj cc j scalar of the metric 

a(t) 2 e 2 ^^hf( X )) by 

[^(x^a^) 2 ^ 2 ). (5.23) 
The trace part of Einstein's equation's ij component is given by 

- M 2 + 6#C (2) + 3 (C (1) ) 2 + + g [3] ^ (2) ) 

= -2G X X'^ + (<^ x - 2X(G X + G XX X)) ^ + V° * , (5.24) 



where 



V u (t) = [\k x + K XX X - 0G X - 50G XX X - 24>G XXX X 2 



+ 2 (3^° - 8 t \nA Q ) <P(G X + G xx X)^j . (5.25) 
Combining with eq. (5.17), we get 
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where two integration spatial functions have been absorbed into (x) and £(°)(x) respec- 
tively, and £®, E\ and £® are background quantities, defined respectively as 

<j>K x - 2X(G X + GxxX)^j C° + 2G X X(6U°C - C°) - V° (A ) 2 - 3M 2 

'' (/) " " 2M 2 - 2G X XA° ' (5 - 27) 
0K X - 2X(G X + GxxX)) Cl + 2G x X(6HCl - C 3 °) - M 2 

^ = " 4M>-4G X XA0 ' (5 - 28) 

0K X - 2X(G X + GxxX)) C$ + 2G X X(2HC$ - C° 3 ) + \M 2 l 

^ = 4Mi pl -4GxXA° • (5 - 29) 

The traceless part of Einstein's equation's ij component is given by 

A^) + ZHA^) + 3C (1) A«i - ( - i*} Nfl^ = 0, (5.30) 



which gives rise to the solution 



Win ^ 3C7^(x)cW;.(x) f t dt > | MflP)«(x)-^WflW(x) /* 

^ f,xj - «3 y a(t') 3 ^ y 



dt'a(t'), (5.31) 



where again an integration spatial function has been absorbed into j(x) and ^^^^ (x) 

(the Ricci tensor of the metric e 2 ^ (0) ( x )/i^(x)) is related to ^R^ 1 - (the Ricci tensor of the 
metric a(t) 2 e 2 < {0) ^ h\f {*)) by 

®RW(x) = a(t) 2 ®RW i j . (5.32) 

From eq. (3.9), we can derive 



+ 44° ) (x)C«i ( x ) C«J ( x ) /-| 



df 



) 3 J a(t') 3 
t A4II rt" di J 



' I a{t'f 
dt" 



2^ } (x)(^^(x)4^^( 2 )(x)) f-** 



t" 

dt'a(t'), (5.33) 



where an integration spatial function has been absorbed into /^^ (x). The ti component of 
Einstein's equation at the 0(e 2 ) order become constraints for the 0(e) order quantities 

- 2M 2 diC {1) - M%D<MAW = {K x 4> + GHGxX)^ - 2G x Xd i ^ 1 \ (5.34) 

where , of order 0(e) itself, is the covariant derivative associated with the 0(e°) order 
metric 

e 2C(°>(x )/t (0) (x) Thig 

gives rise to 3 constraints on the unspecified integration functions 

c£ } }(*)- 



- 10 - 



5.4 Summary 

Here we summarize the solution obtained up to the 0(e 2 ) order: 



C(*,x) 



<£(*,x) 



A}(i,x) 



2 r* dt" f^dt'S^t') 



+ C«f(x)C^Ux) 
+ [ 3 l J R( 2 )(x) ' <l ' 



a(t") 3 
" di'£ 2 (Oa(O 3 



a(t') 



A3 



a(t") 



dt" 

d(t") 3 J a(t') 6 
'"dt'£§(Oa(*') 3 



a{t'f 



+ 0{e 



(5.35) 



K»w)7 



a(t') 3 

*d^(0 [ 3 l^( 2 )(x) 
a(t') 6 + 2 



a(t') 2 



T(i)i 



a(t') 



/\6 



C^j(x) 3C^(x)C 
a 3 



(1)/ 
C 



+ 0(e 3 

dt' 

a(t') 3 



(5.36) 



+ j±j s j i dt / a(t , ) + 0(e 3 ); 



(5.37) 



^•(t,x) = ^(x)-2^(x)C^(x) 



+ 6^(x)C«(x ) ci 1 )j(x) / 



+ 44 0) (x)ci 1) [(x)4 1) J fc (x) 



f dt' 
a(t') 3 

* dt" 
a(t") 3 
* dt" 



a(t") 3 



24 0) (x) (MflWJ^-^MflW^; 



' dt' 
a(t') 3 

a(t') 3 
* dt" 



a(t 



/A3 



dt'a(t'), 



(5.38) 



where A is defined by eq. (5.6), a(t) is defined by eq. (5.9), C°, Cf, C\ and Cg are defined 
by eqs. (5.18-5.21) respectively, £®, £® and £® are defined by eqs. (5.27-5.29) respectively, 
I 3 li?( 2 )*-(x) and [ 3 li?( 2 )(x) are 3D curvature tensors of the metric e 2 ^ (0> ( x )/i^(x). 

There are several unspecified spatial functions in the general solution: £(°)(x), /ij^(x), 



cf^x), C^ ; (x) and C^}(x). (There are a few other C(e) and C(e ) unspecified spatial 
functions, which have been absorbed into these unspecified functions.) /i-^(x) is symmetric 



t(l), 



and has a unit determinant and C^*-(x) is symmetric and traceless, so they each have 5 

degrees of freedom. 3 components of *-(x) are related to other unspecified spatial functions 
respectively by the constraint equations (5.34). 

To count the physical degrees of freedom we should identify the residual gauge freedom. 
Performing the infinitesimal coordinate transformation 



(5.39) 
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and requiring that the synchronous gauge condition on the lapse N and shift N l are respected, 
one straighforwadly obtains that rf" should be of the form 

if = ?] (x) (5.40) 

rf = y*dtV i (i',x)6i j r7 (x) + rf (x) (5.41) 



(77^ may be chosen as 0(e)). From this we infer that the residual gauge freedom amounts to 

,(o) 



4 functions of space. 3 of those can be chosen so as to eliminate 3 spatial functions in (x) 



and 1 chosen so as to eliminate (x). Therefore, we may count the degrees of freedoms as 
follows: 

£(°)(x) 1 scalar growing mode = 1 component, (5-42) 
h^'ipc) 2 tensor growing modes = 5 components — 3 gauge DoFs, (5.43) 
0£ (x) 1 scalar decaying mode = 1 component, (5.44) 
C^*-(x) 2 tensor decaying modes = 5 components — 3 constraints. (5.45) 



6 Late time of Galileon inflation 

At the end of the last section, we claimed that (x) represents a decaying mode. However, 
this is actually not obvious from the general solution given above. After all, a(t) is not the scale 
factor a(t) but is given by a rather complicated expression in terms of background quantities. 
Additionally, and £® also have time dependence. In this section, we would like to briefly 
re-derive the solution for an important special case, the late time of Galileon inflation. This 
will not only allows us to show explicitly that (x) represents a decaying mode, but it 

will demonstrate how one can eliminate the gauge mode (7^ (x) on inflationary backgrounds. 
Moreover, the assumption of quasi-de Sitter expansion drastically simplifies the solution and 
allows an intuitive understanding of its behaviour. Physically, perturbations coming from the 
late time of inflation are observationally most important, as it is these perturbations that 
seed the large scale structure of the observable Universe. Note that, similar to the previous 
section, we mostly suppress the background quantities' order indication W to simplify the 
equations. 

Eq. (5.3) can be integrated to get 

K x cj) + QHG X X oc a(t)~ 3 , (6.1) 

which is an attractor of the dynamical system. So for the later time of Galileon inflation J* 
essentially vanishes. In this limit, the background equations of motion can be simplified to 

K = -3M^H 2 , (6.2) 
K x = -3G x H<j), (6.3) 

and <f> and H become constant and a oc e Ht [35]. Furthermore, we have 



H° -> H, a(t) -> a(t). (6.4) 



The background quantities defined in the last section now all become constant and can also 
be simplified: 

4>K XX X + 6HG X X + 6HG X xX 2 ' 
C° = C^(A ) 2 + C%A° - 3Cl (6.6) 
c o = -\K X + 4K XX X + 2K XXX X 2 + 21H<PG XX X + 6H<pG X xxX 2 

1 2<])K XX X + 12HG X X + 12HG XX X 2 
go 9G X X + 6G XX X 

2 ~ 4>K X xX + 6HG X X + 6HG XX X 2 ' 

M 2 

C 3 ° = _ * , (6.9) 

2<j)K xx X + 12HG X X + 12#G X x^ 2 

V° = -\k x + K XX X + 6H<j>G xx X, (6.10) 



+ (i°) 2 + 3 M p 2 
2M 2 - 2G X XA° 



^K x C! + M 2 pl 

t2 ~ 4M*-4G X XA°' [ } 

£°3 = \el (6-13) 

Note that we have assumed <j>K xx X +6HG X X +6HG XX X 2 ^ and A +6M 2 H/<j>K x + 0, 
which, by using the background EoMs, is equivalent to Gx{Kx — KxxX) + KxGxxX ^ 
and K(G X K X - 2G x K X xX + 2K x G X xX) + K 2 X G X X ^ 0. So the covariant cubic Galileon 
case is included in our solution. We will not discuss the special cases where any of the 
aforementioned quantities actually vanish, but it is easy to follow our formalism in the last 
section to get the relevant results. The constraint eq. (5.34) now becomes 

2(G X XA° - A#)g,cf (x) = M3D«C<ft(x). (6.14) 
Finally, the solution for the late time of Galileon inflation is given by 

+ 0(e 3 ), (6.15) 



2tf 2 a 2 

0(t,x) = ^°)(t)+c( i) (x)- -y ; +^v 2) (^x) ■ v 



ZHa 3 s v ' 1 6Ha 6 

W ) (x) g^|w^w 

4#a 2 + 12#a 6 + 1 h { ' 

A)(t,x) = + C A + — + 0(e% (6.17) 

h . , 0)( | 2/ff(x)C< 1 >*(x) 2^(x)([3]^) ? (x)-^[3] fl (2) (x) ) 
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| 3^(x)^(x)C^(x) + 2^(x)C^i(x)^(x) 

+ ^6 + (6-18) 

Now, we want to explicitly do away with the gauge mode C9~^ (x) in the case of near de 
Sitter inflation by re-slicing. To this end, we choose 

t = t + r ? °(x), (6.19) 
= x l + T}\ (6.20) 

with 

^(x) 

,°(x) = ^i, (6.21) 

Let us consider the effects of the temporal transformation on <f>(°\t): Taylor expansion 
yields ^°\t) = cp^\t) - C^(x) + C(e 3 ), which straightforwardly removes the constant mode 
in eq. (6.16). Though far less obvious, any other effect of the temporal or the spatial part of 
he transformation leads to corrections that are either 0(e 3 ) or can be absorbed in redefinitions 
of (x) and *-(x). The end result is that by re-slicing one can eliminate (7^ (x) with 
all the other terms in the solution unchanged. 



7 Conclusion 



Non-Gaussianity in primordial perturbations can be a powerful probe of different inflation 
models and the associated fundamental theory on which they are based. To calculate or 
evolve non-Gaussian perturbations, one has to go beyond linear perturbation theory. Com- 
plementary to second order cosmological perturbation theory that is applicable to subhorizon 
perturbations, gradient expansion is often used to tackle superhorizon perturbations. In this 
paper, we have developed the superhorizon gradient expansion formalism for Galileon infla- 
tion, a novel inflation model characterized by its higher order derivative interactions. This 
model is inspired by a new class of infrared modifications of gravity, called (generalized) 
Galileon models, introduced to explain the late time accelerated cosmic expansion. There are 
many interesting features in Galilleon inflation, including new shapes of non-Gaussianty [45]. 
We have solved the equations of motion of Galileon infl ation up to se cond order in gradient 
expansion in the synchronous gauge, and obtained the general solution without imposing ex- 
tra conditions on the first order quantities. We have identified the physical degrees of freedom 
in the solution, taking particular care in keeping track of the residual gauge freedom left after 
imposing the synchronous gauge condition. We have also defined a curvature perturbation 
TIP) conserved up to first order. With the formalism developed here, one can evaluate non- 
Gassianities in Galileon inflation at superhorizon scales, and, combined with the conventional 
non-linear perturbation analysis, one can then use the existing data to constrain the model 
parameters. 

A typical way to generate non-Gaussianity is to have multiple fields. Non-Gaussianity 
in the multi-Galileon model has been discussed [46]. It would be interesting to develop a 
superhorizon gradient expansion for multi-Galileon inflation. Another key feature of Galileon 
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gravity is that it is supposed to give rise to 0(1) corrections to general relativity at large 
distances and yet satisfies stringent constraints at short distances, such as in the solar system 
where any modification is typically constrained below 0(1O~ 5 ). This is achieved due to 
the high degree of non-linearity of the Galileon derivative interactions and the phenomenon 
is called the Vainshtein mechanism, originally discovered in massive gravity [43, 55]. This 
mechanism is not easy to see in perturbation theory due to its non-linear nature, and the full 
non-linear problem is difficult to solve. It would be interesting to use the gradient expansion 
in order to get a deeper understanding of the behaviour of these non-linear interactions, at 
least in the regime where it is applicable. 



A Dependence of equations of motion in general covariant scalar-tensor 
theory 

Consider a general covariant scalar-tensor theory of <j) anci 9uuj given by the action S((f>, g^). 
The equations of motion for this system are 



_ _ _L 6S_ 



£ = —— = 0, £ = =0j (A.l) 



with the variation of the action (modulo boundary terms) is given by 

5S = J dx D ^g-(£5 ( j ) + £^5g^)- (A. 2) 

Now, we assume this action is covariant, which means it is invariant under the following 
transformation (Sgx 1 * = — 

Stf = Ctf = fV^, 5 t gT = CtgT = 2V^" ) - (A.3) 

That is, we have 

= J a D x^g~ {£ ■ ?V V <I> + £^ • 2V^r) = 0. (A.4) 

After integration by parts, we get J dPxy/—g£ > v (£\/ u cj) + 2V' i 5^ !/ ) = 0. Since ^ v is arbitrary, 
we have 

SV v <f> = 2V^. (A.5) 

So if the Einstein equations are satisfied (£nu = 0) , the scalar equation of motion is automat- 
ically satisfied {£ = 0). 

B Gradient expansion of the equations of motion in synchronous gauge 

Here we list the Einstein tensor, the effective energy momentum tensor and the scalar current 
(the t component) up to order 0(e 2 ) in the superhorizon gradient expansion. As in the main 
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text, we suppress the background quantities' order indication (°). For example, </>(°) is written 
as (p 

The quantities at the 0(e°) order: 

4 0) =3^ 2 , (B.l) 

4 0) = 0, (B.2) 

= -5) (2H + 3tf 2 ) , (B.3) 

Tf\ t = -K + 2K X X + 6H<j>G x X, (B.4) 

ifu = 0, (B.5) 

Tf ] ) = K - 2G X X4>, (B.6) 

J*(°) = K x <j) + 6HG X X. (B.7) 

The quantities at the 0(e) order: 

4 1) =6ffC (1) , (B.8) 

g£> = ^Di/C* ) = 0, (B.9) 

G^) = -2 (c (1) + 3HC {1) ) 8) - (SHAW) + A^fj , (B.10) 
Tj 1} « = (j>K x + 2j>K xx X + ISifGxX + 12HG XX X 2 ^ jF* + 6<^G X *C (1) , (B.ll) 

T^u = 0, (B.12) 

T^) = [(J>K X - 2X(G X + G XX X)) j>W - 2G X X'^~\ 5), (B.13) 

J* (1) = (k x + 2K XX X + 6#0(G X + G XX X)^ + 6G X ^C {1) - (B.14) 
The quantities at the 0(e 2 ) order: 

G g) = I ^M/jP) + 6 (C (1) ) 2 + 12i7C {2) - AWM (1), fc ) , (B.15) 

Gg^^DiC^-DfcAWJ, (B.16) 
G (2)i = [3] G (2)i _ ^(2)i _ 3ff A ( 2 )i _ 3 ^(i) A (i)< 

J J J J 3 J 

- (V 2 ) + 6#C (2) + 3 (C {1) ) 2 + \A^A^ 8% (B.17) 
T^\ t = {4>K X + 24>K XX X + 18HG X X + 12#GxxX 2 )</>( 2 ) + 60G X ^C (2) 



+ ^-K x +4 J fC xx X+2^ xxx X / +9^G x + 2m(/)Gxx^+6^GxxxX 

+ (18G x ^ + 12Gxx^ 2 )C (1) (1) , (B.18) 
Tf ] ti = {<j>K x + QHG x X)d^ - 2G x Xd i ft 1 \ (B.19) 



-2G X X^ + (<^x - 2X(G X + G xx X)^j ^ - 2j>(G x + G XX X)4> (1 ^^ 



+ (±K X + K XX X - 4>G X - 54>G XX X - 24>G xxx X 2 ^j (^M) 2 



(B.20) 
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J*( 2 ) = dG x XC {2) + (k x + 2K XX X + 6i^G x + 6H<pGxxX^J <j) {2) 

+ + 4>K XXX X + 3FG X + 15i?GxxX + 6HGxxxX 2 ^j ' 

+ (6(^G X + §4>G X xX) C (1) (1) - (B.21) 
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